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Abstract 

We study localization properties of a one-dimensional disordered system 
characterized by a random non-hermitean hamiltonian where both the ran- 
domness and the non-hermiticity arises in the local site-potential; its real part 
being random, and a constant imaginary part implying the presence of either a 
coherent absorption or amplification at each site. While the two-probe trans- 
port properties behave seemingly very differently for the amplifying and the 
absorbing chains, the logarithmic resistance u = ln(l + R4) where i?4 is the 
4-probe resistance gives a unified description of both the cases. It is found 
that the ensemble-averaged < u > increases linearly with length indicating 
exponential growth of resistance. While in contrast to the case of Anderson 
localization (random hermitean matrix), the variance of u could be orders 
of magnitude smaller in the non-hermitean case, the distribution of u still 
remains non-Gaussian even in the large length limit. 
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The study of interference of waves multiply scattered from a system of scatterers 
with non-hermitean hamiltonians has of late become very fashionable. There are 
two classes of problems in this respect; one in which the non-hermit icity is in the 
nonlocal part |], of the hamiltonuan and the other in which the non-hermit icity is 
in the local part (typically in one-body potentials) 0, |, |, §, |7|, |, g, [Hj 0, 0, 
[TJ], [T5|, [HJ. The first category is understood to represent, among other things, the 
physics of vortex lines pinned by columnar defects where the depinning is achieved 

by a sufficiently high transverse magnetic field (represented by an imaginary 
vector potential). In the second category, it is well-known that an imaginary term 
in the local part of the non-hermitean hamiltonian behaves like a source or a sink 
(depending on the sign). We will be concerned here with non-hermiticity in the 
one-body potential only. In a disordered chain with random but real-valued site- 
potentials, almost all the states are exponentially localized and hence an incident 
wave (~ e %kx ) propagating in the positive x-direction is completely backscattered 
due to the well-known localization effects ||17|| . On the other hand, a purely ordered 
chain with fixed absorbing site-potentials (sinks for particles) is intuitively expected 
to lead to exponential decay of the transmittance (forward-scattering), and so it 
does. While similar intuition may lead to the expectation that the transmittance 
would increase indefinitely if each of the fixed imaginary site-potentials is amplifying 
(source of particles), in actuality it does not. In a recent paper || (referred to as I 
from now on), it was shown by the author that the sample with amplifying sites do 
also behave as a perfect reflector in the large length limit, because the perfect lead 
boundary conditions at the two ends of the chain imposes a real energy spectrum for 



the non-hermitean matrix |18| . In this work, we generalize over our work in I and 
consider the effects of disorder on the coherent amplification/ absorption at each 
site and look for an unified description for both the cases. 

We consider a quantum chain of N lattice points (lattice constant unity), repre- 
sented by the standard single band, tight binding equation: 

(E — t n )c n = V(c n _i + c n+ i). (1) 



This open quantum system is coupled to the external world (two reservoirs at very 
slightly different electrochemical potentials) with two identical semi-infinite perfect 
leads on either side. Here E is the fermionic energy, V is the constant nearest 
neighbor hopping term which is the same in both the leads and the sample, e n is 
the site-energy, and c n is the site amplitude at the nth site. Without any loss of 
generality, we choose e n = in the leads and V = 1 to set the energy scale. Inside the 
sample, we choose e n = e r + irj where the real part e r is random and the imaginary 
part r] is a fixed real number which may be either positive or negative, and i = T- 
The random e r is obtained from an uniform distribution on [-W/2, W/2\. Since for 
an isolated scattering potential with a positive (negative) 77 the wave-vector (k) has 
a positive (negative) imaginary part, the wave (~ e tkx ^ decays (grows) exponentially 
with x. Thus a medium with positive r\ at each site is called an absorbing medium 
and a medium with negative r\ an amplifying medium. The physical reason for such 
a description lies in the fact that the scattering in any real medium is never perfectly 
elastic and that in many cases the deviation from perfectly elastic scattering may 
be described by absorption through other inelastic channels or amplification due to 
enhancement of the wave-amplitude (e.g., population inversion in an active medium) 
of incident particles or waves. The complex transmission amplitude in the ordered 
case was calculated to be 



t A = 




(2) 



where 



c = (e 



■ik 



e ^-l) 2 ) rf=( e -^- i/c e 7 -l) 2 , 



(3) 



and the decay length I/I7 



l a and the wave-vector k s are given by 




(4) 



and 



V = (e 



e 7 )sin/c 



(5) 
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The transmittance or the two-probe conductance T = gi = \tA\ 2 obtained from the 
Eq.(2) is found to decay monotonically (exponentially) towards zero for a set of 
absorbers (77 > 0). But, for a set of amplifiers (77 < 0), #2 increases first to a high 
value but eventually (for large L) decays as t& ~ e~' 7 ' L . Disorder (in e r ) in ID is 
known to give rise to an exponential decay. Thus in the presence of both disorder 
and absorption/ amplification, one expects the conductance to behave for L — > 00 
as 

* ~ t A e~ KL ~ e - (l7l+K)L , (6) 

where 1/k — £<j is the localization length for the disorder (in the real part) problem. 
Thus the asymptotic behavior of the transmittance for both the absorbing and the 
amplifying case is the same (i.e., exponential decay) and the effective localization 
length £ for the disordered complex site-potential case is given by = l^ 1 + £^" 1 . 
This result has been checked here by using a numerical transfer matrix method [[0| 
both for the amplifying and the absorbing case. The differencei in the behavior 
of T between the two cases (sharp decay as opposed to an initial rise with large 
oscillations) appears only upto L ~ £, and disappears for L ^> £. It was shown in 
I for the ordered case with ^ = 00, i.e., £ = l a . This may be explained as follows. 
For concreteness, let us consider an E close to zero (band centre). Then \d\ > \c\. 
Now in the absorbing case 7 > making the exponentially increasing term in the 
denominator dominating from the beginning and hence T decays without much of 
an intereference (and hence oscillation) from the beginning. On the other hand, 
in the amplifying case 7 < 0, and hence the exponentially growing term in the 
denominator cannot dominate until some large enough length scale. Also there is a 
lot of interference between the growing and the decaying terms in the denominator 
in this case resulting into an initial rise in T with strong oscillations upto a crossover 
length L c = L C (E, 77) determined by the approximate equality of the two terms in 
the denominator of Eq.(2), i.e., by 2|7|L C ~ ln|d/c|. For the amplifying chain in the 
presence of disorder (W > 0), another competition to the amplification peak {T max ) 
appears because of localization effects and L C (E, W,rj) ~ 1/2 £ ln|d/c| — > as £ — » 



0, i.e., either one or both of W and \rj\ — > oo. 

It may be noted though that the reflectance or the two-probe resistance R — \r\ 2 , 
where r is the complex reflection amplitude (due to backscattering) behaves differ- 
ently from the transmittance. The asymptotically constant reflectance is reduced 
(Roo < 1) for r] > (absorbing chain), and is amplified (R^ > 1) for 77 < (am- 
plifying chain) by factors depending on the disorder strength (W), \rj\ and E. The 
behavior (on an average) is very similar to that for the ordered case (see the figures 
in I). For a fixed E and rj < 0, disorder not only reduces the peak position (L c ) as 
discussed above, but also reduces the T max as well as the asymptotically amplified 
value of Roo. In Fig.l, we show this for the case of E =0.1, 77 = — 0.1 (the same as 
in the ordered case of I), and for a disorder strength W =2.0. One can check in this 
case that the l a ~ 20, and ^ ~ 25. Thus the effective localization length should be 
£ = (l/l a + l/£d) _1 — 11, and this value matches that obtained from Fig.l. Further 
it will be noted that the peak T max =1.3 appears at L c =12 (compared to the T max ~ 
2800 and L c = 68 ~ 3l a in the ordered case), and the amplified backscattering -Roo= 
6.3 (compared to 1600 in the ordered case). 

The issue of the phase distribution of the reflection amplitude r is quite important 
[§, [L6|, |2D| because the evolution of the phase does have important bearings on the 



evolution of the resistance, these two quantities being coupled to each other. In the 
ordered case (W = 0), the phase distribution is a 5-function at an angle dependent 
on E, 7] (both magnitude and sign) and the length L for a small system. It is 
independent of L and the sign of rj for L^> l a (see the expression for r in I). On the 
other hand for a hermitean disorder case (rj = 0), the phase (0) distribution Pl{4>) 
is typically quite non-uniform. The distribution is uniform only in the special case 
of weak localization (L ~ It evolves with length and approaches the stationary 
limit for L 3> This stationary distribution P(4>) has in general two peaks whose 



strengths increase (while the separation decreases) with W [20]. They approach 
a single, narrow peak only in the asymptotic limit of W — > 00. Obviously in the 
presence of both disorder and non-hermiticity, there should be a competition between 
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the broadening effect on Pl(4>) due to a finite disorder and the ^-function narrowing 
effect due to the absorption/ amplification. While it is generally accepted that a 
non-zero r\ suppreses phase fluctuation || [L6], we find below that the suppression 
even in the stationary limit (L ^> £) is only partial in the sense that the phase 
fluctuation does not become a 5- function unless |?7| is very large. Further we find 
that this stationary distribution is independent of the sign of 77 (whether amplifying 
or absorbing chain). In this regard, we first show in Fig. 2 the evolution of Pl{4>) 
towards its stationary distribution P(<p) for the case of E—1.0, W=2.5, and [ 77 [ =0.2 
(i.e., for an amplifying and an absorbing chain). For this choice, l a ~ 10, and ^ ~ 
12 (in lattice units). Thus the competition between the two effects should be quite 
strong. In Figs. 2(a) to 2(d), we show the evolution of Pl{4>) for L = 3, 7, 20 and 200 
respectively. In all cases the histograms were obtained from 10,000 configurations. 
In the Fig. 2(d), L = 200 3> l a or £<j, and hence all the distributions have become 
stationary (no further change), and as one may note there is no difference between 
the distributions for 77= +0.2 or —0.2 in this limit. Next in Figs. 3(a)-3(d), we 
show the stationary distribution P(<p) as |?7| is increased from 0.5 to 3.0. The quite 
broad stationary distribution in Fig.2(d) for 77 = 0.2 gets continually shrunk at larger 
values of \rj\. But even for a large \rj\ = 3.0 in Fig.3(d) where the decay length scale 
l a ~ 0.83 (less than a lattice unit), the stationary distribution P(4>) is far from a 
5-function. In contrast to a recent work |I| , the phase distribution is nowhere even 
approximately Gaussian (particularly because of two peaks in general). Further our 
results differ from another recent work || in that P{4>) does not break apart into 
two narrow peaks as \t]\ — > 00. 

Since the two-probe properties behave differently for absorbing and amplifying 
chains and our main purpose here is to look for an unified description, we study the 
evolution of the four-probe resistance R4 = R/T. Another serious reason for doing 
this is that the probability distribution of T or R for the amplifying chains is not 
well-behaved in the sense that all the moments of T or R (including the first one, 



e.g., < T >) diverge even for a finite L [|1^, [L3| . Since disorder plays an important 
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role, we look at the logarithmic resistance in the form of u(L) = ln(l + R4). We 
note that if Ab is the absorption coefficient then R + T + Ab = 1, and if Am is the 
amplification coefficient then R + T — 1 = Am. Thus while the quantity u is nothing 
but the negative of the logarithmic transmittance — InT in the hermitean disorder 
case, u = ln(l — Ab) — InT in the absorbing case and u = ln(l + Am) — InT in the 
amplifying case. For our calculations, we choose the same parameters as in Figs. 2 
and 3, and show in Fig. 4 the evolutions of u(L) and var(u) =< u 2 > — < u > 2 
in double-logarithmic plots. In the Fig. 4 (a) for a hermitean case (r]= 0), both 
the < u > and var(u) diverges algebraically but with different exponents in the 
regime < L < L2 (where L2 ~ £,d) in accordance with a recent proposal by the 
author [Ell regarding the existence of a two-parameter scaling in that regime and 
that in the strong localization limit L ^> they diverge with a single exponent 
(unity) in accordance with the one-parameter scaling theory. In Figs.4(b)-4(d), we 
introduce increasing amounts of non-hermiticity 77 = ±(0.2 - 3.0) and find some 
profound changes in the relative behavior of the moments. The first thing to note 
is that the moments of u evolve with L almost identically for both the absorbing 
and the amplifying cases; the difference cannot be shown in the scale of the figure. 
This already gives the hint that u is the right variable to look at for an unified 
description of both of these cases. In contrast to the two-probe properties, these 
moments of u do not diverge and there are no odd-even (in L) oscillations. Further 
since there are two length scales involved in this problem, one sees a crossover in 
Fig. 4(b) for |?7| = 0.2 from absorption/ amplification dominated (for L < l a = 10 in 
this example) growth in R4 to the growth affected by localization effects as well (for 
L > £d = 12). For larger \rj\ this crossover is still there at smaller lengths and may 
not be clearly discernible. Further, in contrast to the hermitean disorder case where 
var(u) crosses < u > from the lower side (mildly fluctuating, weakly localized) to the 
higher side (fluctuations dominated, strongly localized) in the thermodynamic limit, 
the asymptotic var(u) in the non- hermitean case is always less than < u >. Thus 
the fluctuations in u are getting suppressed by coherent amplification/ absorption 
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and this may be interpreted as a destruction of the effect of localization on var(u) 
(but not on < u >). Indeed for a large \rj\, var{u) may be orders of magnitude 
smaller than < u >. An example is the case of Fig. 4(d) for \rj\— 3.0, where var(u) 
is about 2.5 orders of magnitude smaller than < u >. It clearly demonstrates how 
strongly the coherent amplification/ absorption works against the localization effects 
in destroying the fluctuations in transport. Yet at the same time, we must notice that 
this destruction is never strong enough to bring back the self-averaging character 
(or, diffusive behavior) to the resistance R4. Since this destruction mechanism is an 
altogether coherent effect, it cannot destroy the the exponential growth of the mean 
and the variance as an incoherent effect (say, due to phonons) would have done. 

To check if Pl{u) is at least asymptotically normal (i.e., P(Ra) log-normal), one 
may calculate the kurtosis parameter K = 1/2(3 — c^/c^) where C2 = var{u) and 
C4 = (u— < u >) 4 at each L. For a normal distribution, C4 = 3c| and hence K = 
0. In Fig.5, we plot K(L) as a function of the system size L for the parameters of 
Fig. 4(b). It may not be surprising that K(L) is not close to zero (Pl(u) far from 
normal) for small lengths, but for much larger lengths L ^> 10£ (where £ ~ 5), 
K(L) still fluctuates within [+0.1,-0.1]. Thus u does not seem to be asymptotically 
normal, and hence i? 4 does not seem to be asymptotically log-normal. 

To summarize, we have looked at the effect of backscattering in a disordered 
chain with coherently amplifying or absorbing site-potentials. Disorder enhances 
the exponential decay of transmittance in the presence of absorbers, while it sup- 
presses the transmittance peak (and the amplified asymptotic reflectance) due to 
the amplifiers. The stationary (L — > 00) probability density of the phase of the 
reflection amplitude is generally an asymmetric double-peaked function (not even 
approximately Gaussian) whose width is progressively suppressed towards zero (5- 
function-like form as in the case of a pure amplifying/ absorbing chain) by larger 
amplifying/ absorbing strength \rj\ or by larger disorder. The four-probe resistance 
i? 4 gives an unifying description for both the amplifying and the absorbing chain, 
and unlike in the case of Anderson localization the amplifying/ absorbing effect can 
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suppress the fluctuations of the logarithmic resistance by orders of magnitude com- 
pared to the average logarithmic resistance. Study of the kurtosis indicates that the 
probability density of the resistance -R4 is not asymptotically log-normal. 

The author acknowledges the support and the warm hospitality of the Condensed 
Matter Research Group at the International Centre for Theoretical Physics, Trieste, 
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Figure Captions: 

Fig.l Evolution of the averages of the logarithmic reflectance (<lnR >) and the 
logarithmic transmittance (<lnT >) as a function of chain length L, for a coherently 
amplifying chain with rj = —0.1, Fermi energy E = 0.1, and a disorder strength W 
= 2.0. The averages are calculated at each L for 10000 configurations. 

Fig. 2 Evolution of the distribution of the phase Pl(4>) with L for both an ampli- 
fying and an absorbing chain rj = ±0.2, E = 1.0, W = 2.5. Histograms representing 
the distribution are drawn using 10000 configurations each for (a) L — 3, (b) L — 7, 
(c) L = 20, and (d) L = 200. The distributions for the pure amplifying/ absorbing 
as well as the hermitean (rj = 0) disorder cases are also shown for each L. The 
Fig. 2(d) represents the stationary distribution P(</>). 

Fig. 3 Stationary phase distribution P(4>) for E — 1.0 and W = 2.5 using L = 
200 for (a) \q\ = 0.5, (b) \r]\ = 1.0, (c) \r]\ = 2.0 and (d) \r]\ = 3.0. 

Fig. 4 Evolution of the mean and variance of u = ln(l + i? 4 ), where R4 is the 
four-probe resistance, as a function of L. Both amplifying and absorbing chains 
with 6000 configurations were used, no distinction between them can be made to 
the scale of the figure. The parameters used are E = 1.0, W = 2.5, and (a) \rj\ = 
0.0, (b) M = 0.2, (c) M = 1.0 and (d) \r]\ = 3.0. 

Fig. 5 The dimensionless kurtosis (see text) parameter as a function of length L 
using E = 1.0, W = 2.5, \rj\ = 0.2 and 6000 configurations. 
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Fig.2a: Dist. of phase, L=3, E=1.0, 10000 config. 
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Fig.2b: Dist. of phase, L=7, E=1 .0, 1 0000 config. 
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: Dist. of phase, L=20, E=1.0, 10000 config. 
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Fig.2d: Dist. of phase, L=200, E=1.0, 10000 config 
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Fig. 3a: Stn. Phase Dist., L=200, E=1.0, 10000 config. 
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Fig. 3b: Stn. Phase Dist., L=200, E=1.0, 10000 config. 
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Fig.3c: Stn. Phase Dist., L=200, E=1.0, 10000 config. 
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Fig. 3d: Stn. Phase Dist., L=200, E=1.0, 10000 config. 
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